function res = trial1(omega,x,h)
res1 = FD_cos(omega,x,h);
res2 = FD_cos(omega,x+h,h);
res3 = FD_cos(omega,x+2*h,h);
res = (res3-2*res2+res1)/h^2;
